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ABSTRACT 


We present cosmological constraints from the sample of Type Ia supernovae (SN Ia) discovered and 
measured during the full five years of the Dark Energy Survey (DES) Supernova Program. In contrast 
to most previous cosmological samples, in which supernovae are classified based on their spectra, we 
classify the DES supernovae using a machine learning algorithm applied to their light curves in four 
photometric bands. Spectroscopic redshifts are acquired from a dedicated follow-up survey of the host 
galaxies of the SNe. After accounting for the likelihood of each supernova being a SN Ia, we find 1635 
DES supernovae in the redshift range 0.10 « z « 1.13 that pass quality selection criteria and can 
be used to constrain cosmological parameters. This quintuples the number of high-quality z > 0.5 
supernovae compared to the previous leading compilation of Pantheon+, and results in the tightest 
cosmological constraints achieved by any supernova data set to date. To derive cosmological constraints 
we combine the DES supernova data with a high-quality external low-redshift sample consisting of 194 
SNe Ia spanning 0.025 < z < 0.10. Using supernova data alone and including systematic uncertainties 
we find Qm =0.352 + 0.017 in a flat ACDM model, and (Qm, w) —(0.264*0075.—0.80* 0-15) in a flat 
wCDM model. For a flat wow, CDM model, we find (Qm, wo, Wa) =(0.49573°033,—0.3670-38 —8.8* 12), 
consistent with a constant equation of state parameter to within ~ 2c. Including Planck Cosmic 
Microwave Background data, SDSS Baryon Acoustic Oscillation data, and DES 3 x 2-point data gives 
(Qm, w) =(0.321 + 0.007,—0.941 + 0.026). In all cases dark energy is consistent with a cosmological 
constant to within approximately 2c. In our analysis, systematic errors on cosmological parameters are 
subdominant compared to statistical errors; these results thus pave the way for future photometrically 
classified supernova analyses such as those planned for the Vera C. Rubin Observatory’s Legacy Survey 


of Space and Time. 


Keywords: supernovae, cosmology, dark energy 


1. INTRODUCTION 


The standard cosmological model posits that the en- 
ergy density of the Universe is dominated by dark com- 
ponents that have not been detected in terrestrial exper- 
iments and thus do not appear in the standard model of 
particle physics. Known as cold dark matter and dark 
energy, their study represents an opportunity to deepen 
our understanding of fundamental physics. 

The Dark Energy Survey (DES) was conceived to 
characterize the properties of dark matter and dark en- 
ergy with unprecedented precision and accuracy through 
four primary observational probes (The Dark En- 
ergy Survey Collaboration 2005; Bernstein et al. 2012; 
Dark Energy Survey Collaboration 2016; Lahav et al. 
2020). One of these four probes is the Hubble dia- 
gram (redshift-distance relation) for Type Ia supernovae 
(SNe Ia), which act as standardizable candles (Rust 
1974; Pskovskii 1977; Phillips et al. 1999) to constrain 
the history of the cosmic expansion rate. To imple- 
ment this probe, the DES SN survey was designed to 
provide the largest, most homogeneous sample of high- 


redshift supernovae ever discovered. The two papers 
that first presented evidence for the accelerated expan- 
sion of the universe (Riess et al. 1998; Perlmutter et al. 
1999) used a total of 52 high-redshift supernovae with 
sparsely sampled light-curve measurements in one or 
two optical passbands. Building on two decades of sub- 
sequent improvements in SN surveys and analysis, we 
present here the cosmological constraints using the full 
5-year DES SN dataset, consisting of well-sampled, pre- 
cisely calibrated light curves for 1635 new high-redshift 
supernovae observed in four bands g, r, i, z. 

For the last decade, SN Ia cosmology constraints 
have largely come from combining data from many sur- 
veys. The recent Pantheon+ analysis (Scolnic et al. 
2022; Brout et al. 2022a) combined three separate mid- 
z samples (0.1 « z « 1.0), 11 different low-z samples 
(z « 0.1), and four separate high-z samples (z > 1.0), 
each with different photometric systems and selection 
functions (Gilliland et al. 1999; Hicken et al. 2009; Riess 
et al. 2001, 2004, 2007; Sullivan et al. 2011; Hicken et al. 
2012; Suzuki et al. 2012; Ganeshalingam et al. 2013; Be- 
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toule et al. 2014; Krisciunas et al. 2017; Foley et al. 2017; 
Riess et al. 2018; Sako et al. 2018; Brout et al. 2019b; 
Smith et al. 2020a). The DES sample, which rivals in 
number the entirety of Pantheon+, does not have the 
low-redshift (z < 0.1) coverage to completely remove 
the need for external low-z samples, but at higher red- 
shift enables us to replace a heterogeneous mix of sam- 
ples with a homogeneous sample of high quality, well- 
calibrated light curves. 

A key aim of the DES analysis was to minimize sys- 
tematic (relative to statistical) errors to enable a robust 
analysis. Vincenzi et al. (2024) shows that our error 
budget is dominated by statistical uncertainty, in con- 
trast to most SN cosmology analyses of the last decade, 
for which the systematic uncertainties equalled or ex- 
ceeded the statistical uncertainties (Betoule et al. 2014; 
Scolnic et al. 2018; Dark Energy Survey Collaboration 
2019). We also highlight that the most critical sources 
of systematics are those related to the lack of a homo- 
geneous and well calibrated low-z sample. 

As the DES sample enables a SN Ia measurement of 
cosmological parameters that is largely independent of 
previous SN cosmology analyses, we have been careful 
to “blind” our analysis. The analysis work described in 
Vincenzi et al. (2024), which stops just short of con- 
straining cosmological parameters, was shared widely 
with the DES collaboration, evaluated, and approved 
before unblinding. Unblinding standards included mul- 
tiple validation checks with simulations and full account- 
ing and explanation of the error budget. No elements of 
the analysis were changed after unblinding. 

In this paper we review the analysis of the complete 
DES SN dataset (as detailed in many supporting pa- 
pers; see Fig. 1) and present the cosmological results. 
An important advance on most previous analyses is that 
we use a photometrically classified rather than spectro- 
scopically classified sample (Moller & de Boissiére 2020; 
Qu et al. 2021), and implement advanced techniques to 
classify SN Ia and incorporate classification probabilities 
in the cosmological parameter estimation (Kunz et al. 
2012; Hlozek et al. 2012). While this advance increases 
the complexity of the analysis, in this work and previ- 
ous papers (Vincenzi et al. 2023; Moller et al. 2022) we 
show that the impact of non-SN Ia contamination due 
to photometric misclassification is well below the statis- 
tical uncertainty on cosmological parameters, and this 
constitutes one of the key results of our analysis. 

Combining our DES data with a low-redshift sample 
(see Sect. 2), we fit the Hubble diagram to test the stan- 
dard cosmological model as well as multiple common ex- 
tensions including spatial curvature, non-vacuum dark 
energy, and dark energy with an evolving equation of 


DES-SN5YR analysis overview 


Data: 
- Calibration (Burke et al. 2018, Brout et al. 2022, Rykoff et al. 2023) 
- SN photometry (Brout et al. 2019, Sanchez et al. 2024) 

- SN spectroscopy (Smith et al. 2020a) 

- DCR and chrom (Lasker et al. 2018, Lee& Acevedo et al. 2023) 

- Host galaxy redshifts and properties (Lidman et al. 2020, Carr et 
al. 2021, Wiseman et al. 2020/2021, Kelsey et al. 2023) 


Simulations: 
- Survey selection effects (Kessler et al. 2019a, Vincenzi et al. 2020) 

- SN Ia intrinsic and dust properties (Brout&Scolnic 2021, Popovic 
et al. 2021a/b, Wiseman et al. 2022) and rates (Wiseman et al. 2021) 
- Contamination (Vincenzi et al. 2019/2020, Kessler et al. 2019b) 


Analysis: 

Pipeline and Overview (Hinton et al. 2020, Vincenzi et al. 2024) 

- Light-curve fitting (Taylor et al. 2023) 

-SN classification (Möller & de Boissière 2020, Qu et al. 2021, 
Vincenzi et al. 2021, Moller et al. 2022) 

- “BEAMS” and bias corrections (Kessler & Scolnic 2017), unbinning 

the SN Hubble diagram (Brout et al. 2020, Kessler et al. 2023) 

- Effects of host galaxy mismatch (Qu et al. 2023) 

- Cosmological contour validation (Armstrong et al. 2023) 


Cosmological results: DES Collaboration 2024 


Testing non-standard cosmological models (Camilleri et al. 2024) 


Figure 1. Overview of supporting papers for DES-SN5YR 
cosmological results. 


state parameter. In Camilleri et al. (in prep. 2024) we 
present fits to more exotic models. 

The structure of the paper is as follows. We begin in 
Sec. 2 by describing the dataset, its acquisition, reduc- 
tion, calibration, and light-curve fitting. We summarize 
the models we test in Sec. 3 before presenting the results 
in Sec. 4; our discussion and conclusions follow in Sec. 5 
and Sec. 6. The details of our data release, which in- 
cludes the code needed to reproduce our results, appear 
in Sánchez (in prep. 2024). 


2. DATA AND ANALYSIS 
2.1. DES and Low-redshift SNe 


Our primary dataset is the full five years of DES SNe, 
which we combine with a historical set of nearby su- 
pernovae from CfA3 (Hicken et al. 2009), CfA4 (Hicken 
et al. 2012), CSP (Krisciunas et al. 2017, DR3) and the 
Foundation SN sample (Foley et al. 2017). "We refer 
to the combined DES plus historical dataset as DES- 
SN5YR. 

'The DES supernova program was carried out over five 
seasons, August to February from 2013-2018, during 
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Figure 2. All DES light curves, showing observed magnitudes in g, r, i, and z bands (left to right respectively) normalized 
by the maximum brightness of each light curve, and with the time-axis de-redshifted to the rest-frame. Each light curve has 
been arbitrarily offset by their redshift, with higher-redshift objects higher on the plot (as labeled on vertical axis). Lines show 
best-fit SALT3 light-curve fits. The g-band and r-band light curves are not used above z ~ 0.4 and z ~ 0.85 respectively because 
that corresponds to the redshifts at which the lower-wavelength limit of the SALT3 model (3500A in the rest frame) passes out 


of their observed wavelength ranges. 


which we observed ten ~ 3deg? fields with approxi- 
mately weekly cadence in four bands (g,r,i,z). Eight 
of the fields were observed to 5e depth of ~23.5 mag in 
all four bands (shallow fields) and two to a deeper limit 
of ~24.5 mag (deep fields). See Flaugher et al. (2015) 
for a summary of the Dark Energy Camera, Smith et al. 
(2020a) for a summary of the supernova program, and 
Diehl et al. (2016, 2018) for observational details. 

'The DES SNe were discovered via difference imaging 
(Kessler et al. 2015) based on the method of Alard & 
Lupton (1998). DES images are calibrated following 
the Forward Global Calibration Method (FGCM; Burke 
et al. 2018; Sevilla-Noarbe et al. 2021; Rykoff 2023), and 


both DES and low-z samples are recalibrated as part 
of the SuperCal-Fragilistic cross calibration effort de- 
scribed in Brout et al. (2022b). SN fluxes are determined 
using scene modeling photometry (Brout et al. 2019b); 
we include corrections from spectral energy distribution 
variations (Burke et al. 2018; Lasker et al. 2019) and 
from differential chromatic refraction and wavelength- 
dependent seeing (Lee & Acevedo et al. 2023). We esti- 
mate the overall accuracy of our calibrated photometry 
to be S 5 mmag. Host galaxies are assigned follow- 
ing the directional light radius (DLR) method (Sullivan 
et al. 2006; Gupta et al. 2016; Qu et al. 2023), and host 
galaxy properties are determined as described by Kelsey 
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Figure 3. Histogram showing the redshift distribution of 
the DES-SN5YR sample, with new DES SNe in blue and 
our low-z sample in red. For comparison the distribution of 
redshifts in the existing Pantheon4- sample is shown in grey 
(Brout et al. 2022a), which also includes the DES SNe from 
the DES-SN3YR analysis (blue dashed line). The five-year 
DES sample contains ~ 4x more supernovae above z ~ 0.4 
than the Pantheon-- compilation. 


et al. (2023) based on Fioc & Rocca-Volmerange (1999) 
using deep coadded images by Wiseman et al. (2020). 
Host galaxy spectroscopic redshifts are obtained primar- 
ily within the OzDES programme (Yuan et al. 2015; 
Childress et al. 2017; Lidman et al. 2020). The final 
data release of photometry of ~ 20, 000 candidates, red- 
shifts of hosts, and host galaxy properties is presented 
in Sánchez (in prep. 2024). 

We apply strict quality cuts to this sample of can- 
didates to select our final high-quality sample for the 
Hubble diagram. The same quality cuts were applied to 
both the low-z sample and the DES supernovae. First, 
we require a spectroscopic redshift of the host galaxy, 
good light-curve coverage (at least two detections with 
SNR» 5 in two different bands), and a well converged 
light-curve fit using the SALT3 model (Kenworthy et al. 
2021; Taylor et al. 2023); this reduces the DES sam- 
ple size to 3621. Additional requirements include light- 
curve parameters (stretch and colour) within normal 
range for SNe Ia, a well constrained time of peak bright- 
ness, good SALT3 fit-probability, and valid distance-bias 
correction from our simulation (see Table 4 of Vincenzi 
et al. 2024, for more detail). Our final Hubble-diagram 
sample includes 1635 supernovae, of which 1499 have a 
machine-learning probability of being a Type Ia greater 
than 5096 (see Sec. 2.2). Note that we do not perform a 
cut on this machine-learning probability, rather we use 
it in the BEAMS formalism that produces our Hubble 
diagram and to weight the SN distance uncertainties in 
the fits to the final Hubble diagram (Kessler et al. 2023). 
The set of all DES light curves is visualised in Fig. 2. 

Since we focus on minimizing potential systematic er- 
rors, we only use the best-calibrated, most homogeneous 
sample of low-z SNe Ia. To reduce the impact of peculiar 


velocity uncertainties we remove SNe with z « 0.025. 
We furthermore combine only a subset of the available 
low-redshift samples: CfA3&4, CSP, and Foundation 
SNe, which are the four largest low-z samples with the 
most well-understood photometric calibration. Our low- 
z sample thus totals 194 SNe with z « 0.1; this can be 
compared to Pantheon+, for which the low-z sample 
was almost four times larger (741 SNe at z « 0.1). We 
have thus exchanged the statistical constraining power 
of more low-z SNe for better control of systematics. The 
redshift distribution of our sample compared to the com- 
pilation of historical samples in Pantheon+ is shown in 
Fig. 3. To conclude, the final DES-SN5YR sample in- 
cludes 1635 DES SNe and 194 low-z external SNe, for a 
total of 1829 SNe. 


2.2. From light curves to Hubble diagram 


A critical step in the cosmology analysis is to con- 
vert each supernova's light curve (magnitude vs time in 
multiple bands; see examples in Fig. 2) to a single cali- 
brated number representing its standardized magnitude 
and estimated distance modulus. 

To achieve this, we use the SALT3 light-curve fitting 
model as presented in Kenworthy et al. (2021); Taylor 
et al. (2023) and retrained in Vincenzi et al. (2024) to de- 
termine the light-curve fit parameters, amplitude of the 
SN flux (xo), stretch (xı), and color (c). These fitted 
parameters are used to estimate the distance modulus, 
1L zm — M, using an adaptation of the Tripp equation 
(Tripp 1998) that includes a correction for observed cor- 
relations between SN Ia luminosity and host properties, 
YGhost- This correction has historically been described 
as a “mass step" but we also consider the possibility 
that it is a “color step" (see Sec. 2.2 of Vincenzi et al. 
2024), 


Hobs,i = mai d- O1, — BC: 4- Y Ghost; — M — Appias,i, (1) 


where m, = —2.5log;o(zo).! The constants a, 8, and 
y are global parameters determined from the likelihood 
analysis of all the SNe on the Hubble diagram, while the 
terms subscripted by i refer to parameters of individual 
SNe. We find a = 0.161 + 0.001, 6 = 3.12 + 0.03, and 
^ = 0.038 + 0.007. We marginalize over the absolute 
magnitude M (see Sect. 3). The final term in Eq. 1 
accounts for selection effects, Malmquist bias, and light 
curve fitting bias. 

The nuisance parameters and Appias,; term in Eq. 1 
are determined using the “BEAMS with Bias Correc- 


1 Following Marriner et al. (2011), we replace the traditional mp 
notation with mz, because in the SALT2 and SALT3 models the 
amplitude term, zo, is not related to any particular filter band. 
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Figure 4. Hubble diagram of DES-SN5YR. We show both the single SN events and the redshift-binned SN distance moduli. 
Redshift bins are adjusted so that each bin has the same number of SNe (~ 50). The 1635 new DES supernovae are in blue, 
and in the upper panel they are shaded by their probability of being a Type Ia; most outliers are likely contaminants (pale 
blue). The inset shows the number of SNe as a function of redshift (same z-range as the main plot). The lower panel shows 
the difference between the data and the best fit Flat-wCDM model from DES-SN5YR alone (third result in Table 2), and 
overplots three other best fit cosmological models — Flat-ACDM model from DES-SN5YR alone (magenta line, first result in 
Table 2), Flat-wowaCDM model from DES-SN5YR alone (green line, fourth result in Table 2), and Planck 2020 Flat-ACDM 


model without SN data (dashed line, QE/^"** =0.317 + 0.008). 


tions” (BBC) framework (Kessler & Scolnic 2017). In 
particular, bias corrections Ajipias,; are estimated from 
a large simulation of our sample. The simulation mod- 
els the rest-frame SN Ia spectral energy distribution 
(SED) at all phases, SN correlations with host-galaxy 
properties, SED reddening through an expanding uni- 
verse, broadband griz fluxes, and instrumental noise 
(see Fig. 1 in Kessler et al. 2019a). Using Eq. 1 there re- 
mains intrinsic scatter of ~ 0.1 mag in Hubble residuals. 
Following the numerous recent studies on understanding 
and modelling SN Ia dust extinction and progenitors 
(Wiseman et al. 2021, 2022; Duarte et al. 2022; Dixon 
et al. 2022; Chen et al. 2022; Meldorf et al. 2023), we 
model this residual scatter using the dust-based model 
from Brout & Scolnic (2021); Popovic et al. (2023a), 
which improves on the previous commonly used models 
in Kessler et al. (2013) that are based on SALT2 error 
models in Guy et al. (2010); Chotard et al. (2011). This 
intrinsic scatter remains the largest source of systematic 
uncertainty from the simulation. 


As we do not spectroscopically classify the SNe and 
thus expect contamination from core-collapse (CC) su- 
pernovae, we perform machine learning light-curve clas- 
sification on the sample following Vincenzi et al. (2023); 
Moller et al. (2022). We implement two advanced ma- 
chine learning classifiers, SuperNNova (Möller & de 
Boissière 2020) and SCONE (Qu et al. 2021) and use 
state-of-the-art simulations to model contamination (es- 
timated to be ~ 6.5%, see Table 10 and Sect. 7.1.5 of 
Vincenzi et al. 2024). Classifiers are trained using core- 
collapse and peculiar SN Ia simulations based on Vin- 
cenzi et al. (2021) and using state-of-the-art SED tem- 
plates by Vincenzi et al. (2019); Kessler et al. (2019b). 
These DES simulations are the first to robustly repro- 
duce the contamination observed in the Hubble residuals 
(Vincenzi et al. 2021; Vincenzi et al. 2024, Table 10). 

For each SN, the trained classifiers assign a probability 
of being a Type Ia, and these probabilities are included 
within the BEAMS framework to marginalize over core- 
collapse contamination and produce the final Hubble Di- 
agram (Kunz et al. 2012; Hlozek et al. 2012). The final 
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DES-SN5YR Hubble diagram is shown in Fig. 4 and 
includes 1829 SNe. 

As discussed in Kessler et al. (2023); Vincenzi et al. 
(2024), the probability that each supernova is a Type Ia 
(Pia) is incorporated in the BBC fit and used to calcu- 
late a BEAMS probability, Ppitay (see Eq. 9 in Kessler 
et al. 2023). BEAMS probabilities are used to inflate 
distance uncertainties of likely contaminants by a fac- 
tor x 1/,/Pgqa) (see Eq. 10 in Vincenzi et al. 2024). 
Therefore, the released Hubble diagram data in- 
cludes distance bias corrections and inflated dis- 
tance uncertainties (see App. A), enabling users to 
fit the Hubble diagram without applying additional cor- 
rections. With this BEAMS uncertainty weight, we find 
75 SNe with distance modulus uncertainties 9,,i.&nal > 
1 mag and 1331 SNe with c,,;.anaj < 0.2 mag.” 

Vincenzi et al. (2024) stops short of performing cosmo- 
logical constraints but provides the corrected distance 
moduli u along with their uncertainties c,,, redshifts for 
each SN, and a statistical+systematic covariance matrix 
C, which we describe further in Sec. 3. 

Armstrong et al. (2023) presents validation of the cos- 
mological contours produced by our pipeline. Validation 
that our analysis pipeline is insensitive to the cosmo- 
logical model assumed in our bias correction simulation 
appears in Camilleri et al. (in prep. 2024), who also test 
more exotic cosmological models. 


2.3. Unblinding criteria 


Throughout our analysis, cosmological parameters es- 
timated from real data were blinded. We validate our en- 
tire pipeline on detailed catalogue-level simulations and 
examine the cosmological parameters estimated from 
simulations to test that the input cosmology is recov- 
ered. In addition to the many tests described in Vin- 
cenzi et al. (2024), the final unblinding criteria that our 
data passed were: 


e Accuracy of simulations: Reduced x? between 
the distribution of data and simulations across a 
variety of observables (redshift, SALT3 parame- 
ters and goodness of the fit, maximum signal-to- 
noise ratio at peak, host stellar mass) is required 
to be between 0.7 and 3.0 (see Vincenzi et al. 2024, 
Fig. 3-4). 


? Applying a binary classification-based cut (SN Ia or not) is not 
optimal, as it assumes the classification is perfect. However, we 
test the binary-cut-based approach by using only the 1499 SNe 
classified with Pj, > 0.5 and assuming they are a pure SN Ia 
sample, and we show that the measured shift in w is small com- 
pared to the statistical uncertainties (Table 11 of Vincenzi et al. 
2024). 


e Pipeline validation using DES simulations: 
Demonstrate that our pipeline recovers the in- 
put cosmology. We produce 25 data-size simu- 
lated samples (statistically independent) assum- 
ing a Flat-ACDM universe with best-fit Planck 
value of Qm and analyze them the same way as 
real data. We fit each Hubble diagram assuming 
a Flat-wCDM model with a Planck prior and find 
a mean bias of w — Wtrue ~ 0.001 + 0.020, where w 
is the mean value of the marginalized posterior of 
the dark energy equation of state parameter over 
the 25 samples, and wirue = —1 is the model value 
of that parameter input to the simulation. 


e Validation of contours: ensuring that our un- 
certainty limits accurately represent the likelihood 
of the models (Armstrong et al. 2023). 


e Independence of reference cosmology: ensur- 
ing that our results are sufficiently independent of 
cosmological assumptions that enter our bias cor- 
rection simulations (Camilleri et al. in prep. 2024). 


2.4. Combining SN with other cosmological probes 


We combine the DES-SN5YR cosmological constraints 
with measurements from other complementary cosmo- 
logical probes. In particular, we use: 


e Cosmic Microwave Background (CMB) measure- 
ments of the temperature and polarisation power 
spectra (TT TEEE) presented by the Planck Col- 
laboration (2020). We use the Python imple- 
mentation of Planck's 2015 Plik lite (Prince & 
Dunkley 2019). 


e Weak lensing and galaxy clustering measurements 
from the DES3x2pt year-3 magnitude-limited 
(MagLim) lens sample; 3 x 2-point refers to the 
simultaneous fit of three 2-point correlation func- 
tions, namely galaxy-galaxy, galaxy-lensing, and 
lensing-lensing correlations (Dark Energy Survey 
Collaboration 2022, 2023). 


e Baryon acoustic oscillation (BAO) measurements 
as presented in the extended Baryon Oscillation 
Spectroscopic Survey paper (eBOSS; Dawson et al. 
2016; Alam et al. 2021), which adds the BAO 
results from SDSS-IV (Blanton et al. 2017) to 
earlier SDSS BAO data. Specifically, we use 
“BAO” to refer to the BAO-only measurements 
from the Main Galaxy Sample (Ross et al. 2015), 
BOSS (SDSS-III Alam et al. 2017), eBOSS LRG 
(Bautista et al. 2021), eBOSS ELG (de Mattia 
et al. 2021), eBOSS QSO (Hou et al. 2021), and 
eBOSS Lya (du Mas des Bourboux et al. 2020). 


Cosmological Model | Friedmann Equation: E(z) = H(z)/Ho = 


Fit Parameters O 


Flat-ACDM | [Qu(1+ 2)? + (10— Qu)] ^ ae 


ACDM | [Ou +z 


3 Qa -F (1 


( ) 
Flat-wCDM [oma 2+0 
( I'd 


Oy)(1 


Flat-wowa CDM fx 1-£ 


Oy )(1 


Qm - Ug rp] Qm, Qa 
1/2 
gja] Qum, w 
z) tunkua) e bwanaa] Om, Wo, Wa 


Table 1. Variations on the standard cosmological model that are tested in this paper, their Friedmann Equations, and the free 


parameters in the fit. 


When combining these data we run simultaneous 
MCMC fits of the relevant data vectors. We present 
three combinations: the simplest CMB-dependent com- 
bination CMB+SN, a CMB-independent combination 
BAO--3x2pt--SN, and a combination of them all. 


3. MODELS AND THEORY 


We present cosmological results for the standard cos- 
mological model — flat space with cold dark matter and 
a cosmological constant (Flat-ACDM) — and some ba- 
sic extensions, such as relaxing the assumption of spatial 
flatness (ACDM), allowing for constant equation of state 
parameter (w) of dark energy (Flat-wCDM), and in- 
cluding a linear parameterisation for time-varying dark 
energy (Flat-wowaCDM) in which the equation of state 
parameter is given by w = wo + wo (1— a) (Chevallier & 
Polarski 2001; Linder 2003). 

To calculate the theoretical distance as a function of 
redshift we begin with the comoving distance, 


c [^ dz 


Rox(2) = Hs Jy EG)’ (2) 


where Z is the redshift due to the expansion of the 
Universe, E(z) = H(z)/Ho is the normalized redshift- 
dependent expansion rate and is given for each cos- 
mological model by the expression in Table 1, Ro = 
c/(Ho4/|QOx|) is the scale factor with dimensions of 
distance (where subscript 0 indicates its value at the 
present day), and Qk = 1 — Oy — Qa is the curvature 
term. The dimensionless scale factor (a = R/Ro) at the 
time of emission for an object with cosmological redshift 
Zisa=1/(1+2Z). The luminosity distance is given by, 


Di (zobs, Z) = (1 ar Zobs) Ro Sk (x(2)), (3) 


where Zobs is the observed redshift, and the curvature 
is captured by S(x) = sinx, x, and sinh x for closed 
(Qk < 0), flat (Ok = 0), and open (Qg > 0) universes 
respectively.? 


3 When Qk = 0 the term RoSt(x) becomes Rox and can be cal- 


culated directly from Eq. 2, bypassing the infinite Ro. 


To compare data (Eq. 1) to theory we calculate the 
theoretical distance modulus, which is dependent on the 
set of cosmological parameters we are interested in (O, 
given in the right column of Table 1), 


u(z, 9) = 5log,9(Di(z, ©)/1 Mpc) + 25. (4) 


We compute the difference between data and theory 
for every ith supernova, Aw; = fobs, — u(2;, O), and 
find the minimum of 


x = AC Au. (5) 


where C~! is the inverse covariance matrix (including 
both statistical and systematic errors) of the Ap vector 
(see Sec. 3.6 of Vincenzi et al. 2024). 

'The uncertainty covariance matrix includes a diago- 
nal statistical term (discussed Sec. 2.2) and a systematic 
term. The systematic covariance matrix is built follow- 
ing the approach in Conley et al. (2011) and accounts 
for systematics such as calibration, intrinsic scatter, and 
redshift corrections (see Table 6 of Vincenzi et al. 2024). 
Each element of the covariance matrix expresses the co- 
variance between two of the SNe in the sample. The 
covariance matrix has dimensions of the number of su- 
pernovae Ngne x Ngne and we follow the formalism in- 
troduced by Brout et al. (2021) and Kessler et al. (2023). 

Finally, the absolute magnitude of SNe Ia (M) and 
the Ho parameter (which appears in the luminosity dis- 
tance) are completely degenerate and therefore they 
are combined in the single parameter M = M + 
5logio(c/ Ho). All of our cosmology results are marginal- 
ized over this term. Therefore, the value of Hp has no 
impact on the fitting of our cosmological results, and we 
do not constrain Hg. While M has no impact on cos- 
mology fitting, a precise value is needed to simulate bias 
corrections. The M uncertainty is below 0.01, resulting 
in a negligible impact on bias corrections (Brout et al. 
2022a; Camilleri et al. in prep. 2024). 


4. RESULTS 


With the new DES high-redshift supernova sample 
we can put strong constraints on cosmological models. 
Of particular interest is whether dark energy is consis- 
tent with a cosmological constant or whether its density 
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Figure 5. Constraints on matter density in the Flat- 
ACDM model from DES-SN5YR only (cyan), DES-SN5YR 
combined with CMB constraints from Planck Collabo- 
ration (2020) (blue), and DES-SN5YR combined with 
BAO+3x2pt (orange), and all probes combined (DES- 
SN5YR--BAO--3x2pt and CMB constraints, dark orange). 
CMB constraints only and BAO+3x2pt constraints alone 
are also shown for comparison (dashed and dotted-dashed 
respectively). 


and/or equation of state parameter varies over the wide 
redshift range of our sample. The results of our cosmo- 
logical fits are outlined in this section and summarized 
in Table 2, and their implications are explored in Sec. 5. 

We estimate cosmological constraints using the Cos- 
moSIS framework (Zuntz et al. 2015), the samplers 
emcee for best fits (Foreman-Mackey et al. 2013), and 
PolyChord for tension metrics (Handley et al. 2015),* 
except for fits that include BAO--3x2pt, which are cal- 


4 For each emcee fit we use a number of walkers that is at least twice 
the number of parameters and ensure the number of samples in 
the chain is greater than 50 times the autocorrelation function, T 
(Nsamples/T > 50). For each PolyChord fit, we use a minimum of 
60 live points, 30 repeats, and an evidence tolerance requirement 
of 0.1 (except for ACDM with all datasets combined, for which we 
accepted a slightly weaker tolerance because convergence was too 
slow). When combining with other datasets we run simultanous 
MCMC chains including all relevant data vectors. Flat priors 
that encapsulate at least the 99.796 confidence region were chosen 
in each case. 
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Figure 6. Constraints for ACDM model (non-zero curvature 
allowed) from the DES-SN5YR dataset only (cyan), from 
DES-SN5YR combined with BAO and weak lensing mea- 
surements (orange), and from DES-SN5YR combined with 
CMB measurements (blue). For comparison, we also present 
cosmological constraints from Planck Collaboration (2020) 
only (black dashed). 


culated using PolyChord for both best fit and tensions.” 
For all fits we present the median of the marginalized 
posterior and cumulative 68.27% confidence intervals. 
The chains and code (with the flexibility to test other 
statistical choices) will be made available on github upon 
acceptance of this paper. Figs. 5, 6, 7 and 8 all present 
the joint probability contours for 68.3% and 95.5%. 


4.1. Constraints on Cosmological Parameters 
4.1.1. Flat-ACDM 


For the simplest parameterization, Flat-ACDM, € is 
the only free parameter. We show the probability den- 
sity function (PDF) of this constraint for DES-SN5YR 
in Fig. 5; we measure a value of Qm —0.352 + 0.017. 
We also show the probability distribution of the Planck 
Collaboration (2020) measurement of Offene =0.317 + 
0.008. These are approximately 2c apart, but not in 
significant tension as discussed in Sec 4.2. 


5 The main advantage of emcee is it gives slightly more accurate 
best fit x? than PolyChord. However, we decided the tiny im- 
provement in accuracy was not worth the environmental impact 
(Stevens et al. 2020) of the extra compute time (which was sub- 
stantial for the many-dataset fits). 
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Figure 7. Same as Fig. 6 but for Flat wCDM model. The 
horizontal dotted line marks the equation of state values for 
a cosmological constant, i.e. w = —1. 


Combining DES-SN5YR with Planck CMB gives 
Ow —0.338*0 014, while combining with BAO+3x2pt 
gives QM =0.930") ig: Combining all three gives 
Qm =0.315 + 0.007. Interestingly, the combination of 
all data sets (dark orange in Fig. 5) gives a lower Ow 
than any of the other combinations. The reason can 
be seen in Fig. 6, where all constraints cross the Flat 
Universe line to the upper left of any individual best fit. 


4.1.2. ACDM 


Fitting DES-SN5YR to the ACDM model, we find 
(Qm, 04)—(0.291* 0:065, 0.55 + 0.17), consistent with a 
flat universe (Q4 —0.16 + 0.16); see Fig. 6. Combin- 
ing DES-SN5YR with BAO--3x2pt is also consistent 
with a flat Universe, with uncertainties on Qg reduced 
to ~ +0.034, while the combination with Planck gives 
Oy =0.010 + 0.005. The combination of all three gives 
Or —0.002*0.00. 


4.1.3. Flat-oCDM 


Fitting DES-SN5YR to the Flat-; CDM model, we 
measure (Nm, w) = (0.264* 0075.—0.80* 0-14); see Fig. 7. 
This is consistent with a cosmological constant (within 
2c), although our data favors a w-value that is slightly 
larger than —1. 

The w — Qm contours from SN alone are highly non- 
Gaussian with a curved ‘banana’-shaped degeneracy. 
The best fit value for w or Qy is thus an insufficient 
summary of the SN information, as a small shift along 
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the degeneracy direction can result in large shifts in 
the best-fit values. To address this issue, in Camil- 
leri et al. (in prep. 2024) we introduce a new param- 
eter, Qu(z) = —@/(aH}) = q(H/Ho)?. This combina- 
tion of the deceleration parameter q and the Friedmann 
equation H/Ho follows the curve of the degeneracy in 
the w — Qm plane. Therefore, measuring Qy(z) sum- 
marizes the supernova information in a single, almost 
degeneracy-free value. One has to choose the redshift 
at which one quotes Q(z), to best match the angle of 
the degeneracy for the redshift range of the sample. We 
find Qu (z = 0.2) = —0.340 + 0.032 using DES-SN5YR 
only (see Camilleri et al. in prep. 2024). This Qg value 
can be used to roughly approximate the DES-SN5YR 
results and characterize the constraining power without 
the need for a full fit to the Hubble diagram. 

The degeneracy in the w — Qm plane is bro- 
ken by combining SNe with external probes. 
Combining with Planck, we measure (Qw,w) = 
(0.337*0.012,—0.055*0022), again within 2e of a cos- 
mological constant. Planck alone provides only a loose 
constraint on the equation of state parameter of dark 
energy, wPlenk — .—1.51'9?7: combining with DES- 
SN5YR reduces the uncertainty significantly due to the 
different degeneracy direction, demonstrating the com- 
bined constraining power of these two complementary 
probes. 

Combining DES-SN5YR with BAO+3x2pt we find 
w ——0.922*0035. slightly over 2e from the cosmologi- 
cal constant. This data combination demonstrates that 
these late-universe probes alone provide constraints that 
are consistent with — and of comparable constraining 
power to — the combination of SN and CMB data. The 
full combination of all data sets gives w =—0.941+0.026. 


4.1.4. Flat-wowa CDM 


Fitting DES-SN5YR alone to the Flat-wgw; CDM 
model gives an equation of state that is slightly over 
2c from a cosmological constant, marginally preferring 
a time-varying dark energy (Qm, wo, Wa) —(0.495*0 015, 
—0.36*038. —8.8*7-7): see Fig. 8. 

Combining DES-SN5YR and the CMB, we find 
(QM, Wo, Wa) -—(0.325*0015,—0.73 + 0.11,—1.17* 055), 
which again deviates slightly from the cosmological con- 
stant. 'The same trend is seen when combining with 
BAO+3x2pt and with all data combined. The negative 
Wa means that the dark energy equation of state param- 
eter is increasing with time (sometimes referred to as a 
"thawing" model). 


6 Similar to the Sg parameter used in lensing studies to approxi- 


mate og-Qyy, constraints. 
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Figure 8. Same as Fig. 6 but for wow; CDM model. The dashed crosshairs mark the equation of state values for a cosmological 
constant, i.e. (wo, Wa) = (—1,0). The residuals between the DES-SN5YR best fit Flat-wowaCDM w.r.t. Flat-wCDM model are 


presented in Fig. 4. 


4.2. Goodness of fit and tension 
4.2.1. x? per degree of freedom 


'To assess whether our best fits are good fits we cal- 
culate the x? per degree of freedom for all our dataset 
and model combinations; see the last column of Table 2. 
The x? we use for this test is the maximum likelihood 
of the entire parameter space, not the marginalized best 
fit for each parameter. 

'The number of degrees of freedom is the number of 
data points minus the number of parameters that are 
common to all datasets (i.e., the cosmological param- 
eters of interest). The number of data points added 
by the CMB, BAO, and 3x2pt is respectively 615, 8, 
and 471. Due to our treatment of contamination (by 
inflating the uncertainties of SNe with a low Pia), we 
approximate the effective number of data points in the 
DES-SN5YR sample by >) Paa) = 1735 (rather than 
the total number of data points, 1829). 


Ideally, a good fit should have y?/d.o.f.— 1.0. The 
slightly low x?/d.o.f. for the DES-SN5YR data arises 
because >) Pp(i; only approximates the number of de- 
grees of freedom, and the same behaviour is also seen in 
simulations. 


4.2.2. Suspiciousness 


Suspiciousness, S, (Handley & Lemos 2019) is closely 
related to the Bayes ratio, R,’ and can be used to 
assess whether different datasets are consistent. How- 
ever, while the Bayes ratio has been shown to be prior- 
dependent (Handley & Lemos 2019), with wider prior 
widths boosting the confidence, Suspiciousness is prior 
independent. Therefore, Suspiciousness is ideal for cases 
such as ours where we have chosen deliberately wide 
and uninformative priors (Lemos et al. 2021, Sec. 4.2). 
Jeffreys’ scale (Jeffreys 1961) suggests In $ < —2.5 is 


T Suspiciousness, S, is related to the Bayes ratio R and Bayesian 


information J and is defined as In $ = ln R — ln T. 
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Table 2. Results for four different cosmological models, sorted into sections for different combinations of observational con- 
straints. These are the medians of the marginalized posterior with 68.27% integrated uncertainties (‘cumulative’ option in 
ChainConsumer). For each fit we also show the x? per degree of freedom as a measure of the goodness of fit. 


Ou Ox wo Wa x? /dof 
DES-SN5YR (no external priors) 
Flat-ACDM 0.352 + 0.017 g - - 1649/1734=0.951 
ACDM 0.201*0002 0.16 0.16 : - 1648/1733=0.951 
Flat-wCDM 0.26419:974 à —0.80+9:14 - 1648/1733=0.951 
Flat-woWaCDM .— 0.495*0023 - —0.36* 0.20 —8.8* 2: 1641/1732—0.948 
DES-SN5YR + Planck 2020 
Flat-ACDM TEE re - - : 2237 /2349=0.952 
ACDM 0.359*0 0:5 0.010 + 0.005 : " 2231/2348—0.950 
Flat-wCDM (33003 - (B5 to oss - 2234/2348—0.951 
Flat-ugw,CDM 0.325+9-018 - —0.73+0.11  —117*055 2931/2347=0.951 
DES-SN5YR + SDSS BAO and DES Y3 3x2pt 
Flat-ACDM D330 o cic - - - 2194/2212—0.992 
ACDM 0.327*0012 — 0.030 + 0.034 E : 2194/2211=0.992 
Flat-wCDM UK CERE - —0.922*0-037 - 2188/2211—0.989 


Flat-wows, CDM 0.334 + 0.012 


—0.778*0020 — —0.93*03$  2191/2210—0.992 


DES-SN5YR + Planck 2020 + SDSS BAO and DES Y3 3x2pt 


Flat-ACDM 0.315 + 0.007 g : - 2791 /2828=0.987 
ACDM 03S SDN. O00 re. - - 2825 /2827=0.999 
Flat-wCDM 0.321 + 0.007 g —0.941 + 0.026 - 2785/2827=0.985 


Flat-wowaCDM 0.325 + 0.008 


—0.773*0075  —0.8319:33 2782/2826=0.984 
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Figure 9. Measurements of Suspiciousness (AIn(S)) be- 
tween the DES-SN5YR and Planck 2020 datasets for the 
four models constrained in this paper. Further left indicates 
higher tension where the shaded regions reflect “substantial” 
(yellow) and "strong" (red) evidence of tension according to 
Jeffreys’ scale (Jeffreys 1961). The values and uncertainties 
represent the mean and standard deviation of realizations 
estimating sample variance using the ANESTHETIC software. 


"strong" tension, —2.5 « In S « —1.2 is “substantial” 
tension, and In S > —1.2 indicates the datasets are in 
agreement. 

We determine InS using the ANESTHETIC software 
(Handley 2019), which produces an ensemble of real- 
izations used to estimate sample variance. Results are 


Flat-wow,CDM 


2.5 


-2.5 0 


A(In BE) compared 
to FlatACDM 


Figure 10. Bayesian Evidence difference relative to Flat- 
ACDM (A(In BE)). We present the results for the four dif- 
ferent models tested in this analysis and for the three combi- 
nation of datasets used (DES only in cyan, DES+Planck in 
blue, DES4-BAO--3x2pt in orange). An increase (decrease) 
in A(In B E) indicates that a model is disfavoured (favoured) 
compared to Flat-ACDM. 


quoted using the mean of the ensemble, with the error 
bars reflecting the standard deviation. 

In Fig. 9 we plot the Suspiciousness values for the 
DES-SN5YR data vs Planck 2020 and vs BAO--3x2pt 
data. We find no indication of tension using any of the 
four models investigated in this paper. 
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4.3. Model Selection 


Finally, we use Bayesian Evidence to test whether the 
extra parameters in the more complex models we test are 
warranted, given the data. In Fig. 10, we present the dif- 
ference in the logarithm of the Bayesian Evidence, A(In 
BE), relative to Flat-ACDM for the four different mod- 
els tested in this analysis and for the three combinations 
of datasets used in Fig. 10. 

To evaluate the strength of evidence when comparing 
Flat-ACDM with more complex models, we again use 
Jeffreys’ scale. This empirical scale suggests that A(In 
BE)> 2.5 (and < —2.5) is moderate evidence against 
(in support of) the more complex model, whereas A(In 
BE)> 5 (and < —5) is strong evidence against (in sup- 
port of) the more complex model (for a review of model 
selection in cosmology see Trotta 2008). We note that 
none of the datasets considered in this analysis strongly 
favours cosmological models beyond Flat-ACDM. The 
priors that we choose for model comparison are w € 
(—1.5, —0.5), wa € (—10,10) and Qg € (—0.5,0.5). We 
consider these priors (which determine the penalty for 
more complex models) to be reasonable in terms of gen- 
eral considerations, such as avoiding universes that are 
younger than generally accepted stellar ages (see Section 
5.1.3). We also find the results to be consistent with the 
Akaike Information Criteria, another commonly used 
model comparator. 


5. DISCUSSION 
5.1. The big questions 


5.1.1. Is the expansion of the Universe accelerating? 


Twenty five years ago Riess et al. (1998) found 99.5%- 
99.9% (2.80 to 3.90) evidence for an accelerating Uni- 
verse, by considering the deceleration parameter q — 
(aa)a~? and integrating over the likelihood that qo < 0. 
Importantly they note that since go is measured at the 
present day but the data span a wide range of redshifts, 
qo can only be measured within the context of a model, 
either cosmographic or physically motivated. They used 
the ACDM model, in which qo = Qy/2 — Qa. 

Doing the same with DES-SN5YR data gives 
99.99998% confidence (5.20) that qo < 0 in ACDM, or 
a 2 x 1077 chance that the expansion of the Universe is 
not accelerating. As noted in Section 4.1.3, our confi- 
dence is even higher that the universe was accelerating 
at z ~ 0.2. When we further assume flatness, the con- 
fidence in an accelerating Universe is overwhelming (no 
measurable likelihood for a decelerating Universe) and 
we find qo ——0.530*0017. For more fits of qo using a cos- 
mographic approach see Camilleri et al. (in prep. 2024). 


5.1.2. Is dark energy a cosmological constant? 


As seen in Sec 4.1, a cosmological constant is a good 
fit to our data, but not the best fit. Our best fit 
equation of state parameter is slightly (more than 1c) 
higher than the cosmological constant value of w — —1 
(both for SNe alone and in combination with Planck 
or BAO--3x2pt). Our result agrees with the recent 
result from the UNIONS compilation analyzed with 
the UNITY framework (Rubin et al. 2023) (which ap- 
peared while this paper was under internal review). The 
Pantheon4- result (Brout et al. 2022a) is within lo of 
w = —1, but also on the high side (w = —0.90 + 0.14). 

Furthermore, our analysis slightly prefers a time- 
varying dark energy equation of state parameter when 
we fit for w(a) such that the equation of state parame- 
ter increases with time (again for all data combinations), 
known as a “thawing” model. Model selection, however, 
is inconclusive. 

'The constraints on time-varying w are enabled by the 
wide redshift range of the DES-SN5YR sample. Our 
analysis as described in Vincenzi et al. (2024) gives us 
confidence that systematic uncertainties in this data are 
below the level of our statistical precision. Nevertheless, 
it is important to recognize that (a) the low-z sample is 
the one for which we have the least systematic control 
and (b) the very high-redshift SNe are the ones for which 
bias-corrections are large (7 0.1 mag) and more uncer- 
tain (e.g., accurate estimation of spectroscopic redshift 
efficiency is more challenging as we go to higher red- 
shifts), and for which the uncertainties on the rest-frame 
UV part of the SN Ia spectral energy distribution have 
more impact on SN distances estimations (see also Brout 
et al. 2022a). 

'To test whether our fits are dominated by any partic- 
ular redshift range we ran cosmological fits (a) remov- 
ing low-z data (i.e., DES SNe alone) and (b) removing 
high-z data (i.e., removing ~ 80 SNe at z > 0.85, for 
which we use only two bands; see Fig. 2). Most of the 
cosmological results obtained with the subsamples are 
consistent with the results found for the full sample, 
(based on running the same test on 25 simulated data 
samples). However, we found that removing the low-z 
sample shifts the contours in the Flat-wCDM slightly 
down, which would make the combined fits more con- 
sistent with w = —1. The Flat-wgw4 CDM results are 
stable to sub-sample selection. 

We showed in Vincenzi et al. (2024) that systematic 
uncertainties are sub-dominant to the statistical uncer- 
tainties in our sample. Nevertheless, in the future a new 
low-redshift sample (see Sec. 5.3) would help alleviate 
any remaining doubt about calibration and systematics 
in the existing low-z sample, and an even higher-redshift 


supernova survey would help alleviate any modelling 
concerns by minimizing selection effects even at z ~ 1. 


5.1.3. How old is the Universe? 


One of the issues that the discovery of dark energy 
solved is the age of the Universe (tọ) problem — globu- 
lar cluster age estimates, in combination with high es- 
timates of Ho, were inconsistent with models that were 
not accelerating (VandenBerg et al. 1996; Gratton et al. 
1997; Chaboyer et al. 1998). 

Our results, which favor a dark energy equation of 
state parameter slightly higher than w = —1 would im- 
ply that the age is slightly younger than the age found 
in a Universe where dark energy is a cosmological con- 
stant (for the same values of Ho and present dark energy 
density). 

To calculate the Universe’s age, one needs a value of 
Ho in addition to the best fit cosmological model. Since 
we do not constrain Ho in this analysis, we present our 
measurement of the combination Hoto. In other words, 
we give tg in units of the Hubble time ty = 1/Hg.5 Our 
best-fit DES-SN5YR result in Flat-ACDM would have 
an age of (0.921 —- 0.013)£g. This is ~ 3% younger than 
Planck (t)32"** = (0.950 + 0.007)t 5), corresponding to 
an age difference of approximately —0.4 Gyr. Our best 
fit Flat-wowaCDM model gives an age (0.86 + 0.02)t x, 
about 9% younger than the Flat-ACDM Planck result, 
corresponding to an age difference of approximately 
—1.3 Gyr. Such a young age is unlikely given the age of 
the oldest globular clusters (Valcin et al. 2020; Cimatti 
& Moresco 2023; Ying et al. 2023). In the future, this 
information could be used as a prior to limit the feasible 
range of time-varying dark energy. 


5.1.4. Does our best fit resolve the Hubble tension? 


As pointed out in Planck Collaboration (2020, their 
Sec. 5.4), the only basic extensions to the base Flat- 
ACDM model that resolve the Ho tension are those in 
which the dark energy equation of state is allowed to 
vary away from w = —1. In the wCDM model a phan- 
tom equation of state parameter of w ~ —1.5 would 
help resolve the tension (Di Valentino et al. 2021, their 
Sec. 5.1), and it is clear from Fig. 7 that CMB alone 
actually prefers w « —1. In this model, Planck alone 
does not constrain Hp very tightly, and they refrain 
from quoting a value, (see Table 5 of Planck Collabo- 
ration (2020)), but lower w correlates with higher Ho. 
However, the DES-SN5YR data shows a slight tendency 


8 If Ho = 68 kms-! Mpc-!, tg (68) = 14.38 Gyr. 
If Ho = 73 kms~!Mpc~?, tg (73) = 13.40 Gyr. 
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for w > —1, essentially ruling out this solution within 
wCDM. 


5.2. Comparison with DES-SN3YR and Pantheon+ 


It is informative to compare the results of the pre- 
vious DES-SN3YR analysis (Dark Energy Survey Col- 
laboration 2019; Brout et al. 2019a) with the results 
of the DES-SN5YR analysis presented in this work. 
'The DES-SN3YR analysis included 207 spectroscopically 
confirmed SNe Ia from DES and 127 low-redshift SNe 
from CfA and CSP samples (see also Fig. 3). A fraction 
of those events is in common between both analyses (55 
from low-z external samples and 146 DES SNe).? 

However, the DES-SN3YR analysis differs from the 
analysis presented here in many aspects. The SN Ia 
intrinsic scatter modelling has been significantly im- 
proved (from ‘G10’ and constant Gint floor, to the more 
sophisticated modelling of intrinsic scatter introduced 
by Brout & Scolnic 2021; Popovic et al. 2023a), the 
BBC software has been updated (from BBC ‘5D’ and 
a binned approach, to BBC ‘4D’ and an unbinned ap- 
proach), the z4 — M, correlations have been incorporated 
into simulations (following the work by Smith et al. 
2020b; Popovic et al. 2021), and the light-curve fitting 
model has been updated from the SALT2 model to the 
SALT3 model (see Taylor et al. 2023, for a compari- 
son between SALT2 and SALT3 using the DES-SN3YR 
sample). Finally, the DES-SN3YR analysis did not re- 
quire machine-learning classification and the implemen- 
tation of the BEAMS approach because it is a sam- 
ple of spectroscopically selected SNe Ia. We compare 
the final SN distances in Fig. 11 and find consistent 
results (differences in binned distances are on average 
0.02 mag, even in the redshift ranges where contami- 
nation is expected to be high). The cosmological re- 
sults from DES-SN3YR and DES-SN5YR are consistent 
within uncertainties (when assuming Flat-ACDM, Qm 
are 0.331 -E 0.038 and 0.352-E 0.017 for DES-SN3YR and 
DES-SN5YR respectively, while when assuming Flat- 
wCDM and including CMB priors, w are —0.978 4: 0.059 
and —0.955* 0:022). 


''he other main dataset we can compare to is Pan- 
theon+, which contains a significant amount of indepen- 


9? Not all events included in the DES-SN3YR analysis are included 
in the DES-SN5YR analysis and vice-versa. This is due to the two 
analyses implementing different sample cuts. For example the 
z > 0.025 cut and the requirement for a host-galaxy redshift in 
DES-SN5YR exclude respectively 44 and 29 low-z SNe that were 
in the DES-SN3YR sample. DES-SN5YR also uses a new SALT 
model (which affects the SALT-based cuts), and is restricted to 
SNe that pass selection cuts across all systematic tests (see Table 
4 in Vincenzi et al. 2024). 
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Figure 11. Comparison between Hubble residuals for the 
DES-SN3YR and DES-SN5YR analyses w.r.t. the best fit 
Flat-; CDM for the DES-SN5YR analysis. Hubble residuals 
are binned in redshift and we present the weighted mean 
and standard deviation of the mean in each redshift bin. The 
redshift range covered by the low-z sample is highlighted and 
shown with thick dotted lines. The two DES samples are 
consistent with each other. Note the DES-SN3YR analysis 
only includes spectroscopically confirmed SNe whereas the 
DES sample in the DES-SN5YR analysis consists entirely of 
photometrically identified SNe Ia and extends to higher-z. 


dent data (all the high-z data). The DES sample is on 
average much higher redshift than the Pantheon+ sam- 
ple (see Fig. 3), with over a quarter of the DES-SN5YR 
sample being at high enough redshift (z = 0.64) to probe 
the likely decelerating!? period of the Universe (com- 
pared to 6% in Pantheon+). We show a comparison of 
the contours in Fig. 12. We find very similar constrain- 
ing power between Pantheon+ and DES-SN5YR, and 
the DES-SN5YR value of w is within lo of Pantheon+ 
(Brout et al. 2022a). These analyses are not fully in- 
dependent as a fraction of the low-z sample is shared. 
However, all of the high-z dataset is independent, and 
DES is a photometric sample while Pantheon- is fully 
spectroscopic. The constraints on w are similar between 
DES and Pantheon+ as DES high-z has better precision 
per SN than Pantheon4- and has significantly higher sta- 
tistical power at z > 0.4 (see Fig. 3), but Pantheon 
used 2x more low-redshift SNe (which we do not in- 
clude in order to be able to better control systematic 
uncertainties). 


5.3. DES and Next Generation Supernova Samples 


This analysis has shown that moving from a spectro- 
scopically confirmed sample as done in Dark Energy Sur- 
vey Collaboration (2019) to a photometric sample can 
increase the sample size of well-measured supernovae 
significantly (from 207 DES SNe Ia in DES-SN3YR to 


10 The redshift at which the Universe began accelerating in ACDM 
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Figure 12. Constraints in Flat-wCDM from the DES- 
SN5YR sample, the Pantheon+ sample (with and without 
CMB priors), and the Amalgame sample. The constrain- 
ing power of the DES-SN5YR and Pantheon+ samples is 
comparable and consistent, despite Pantheon+ being a spec- 
troscopic SN Ia sample combining 17 different surveys. The 
*Amalgame" sample includes the SDSS and PS1 photometric 
SN samples (> 1700 intermediate-redshift and high-redshift 
SNe), however it does not include a low-z anchoring sample 
(hence the larger contours). DES-SN5YR and Pantheon+ 
are also combined with CMB constraints (for both we use 
the Planck lite Python implementation presented by Prince 
& Dunkley 2019). The horizontal dotted line marks the equa- 
tion of state values for a cosmological constant. 


> 1600 in DES-5YR), consistent with an analysis of 
Pan-STARRS SNe in Jones et al. (2018). This improve- 
ment arises because photometric classification alleviates 
the bottleneck of limited spectroscopic resources. The 
improvement will increase for future surveys as more 
candidates are discovered, but the available time for 
spectroscopy does not increase commensurately. Impor- 
tantly, the work of Vincenzi et al. (2024) shows that 
systematic uncertainties due to photometric classifica- 
tion are not limiting. Instead, the “conventional” sys- 
tematics of calibration and modeling the intrinsic scatter 
remain the most significant challenges. 

There is potential for further increase of the statistical 
power of the DES sample if one moves to using SNe 
in which a host galaxy spectroscopic redshift was not 
acquired and instead relies on photometric redshifts of 
the SNe and the galaxy. This path was explored by 
Chen et al. (2022) for a subset of DES SNe, namely 


ones that occur in redMaGiC galaxies, and has been 
explored as well for SuperNova Legacy Survey (SNLS, 
Ruhlmann-Kleider et al. 2022) and the Vera C. Rubin 
Observatory Legacy Survey of Space and Time (LSST) 
in Mitra et al. (2023). These analyses show that the use 
of photo-zs do not introduce systematic uncertainties 
to a scale similar to the statistical uncertainties. This 
potential is highlighted by the ~ 2400 SNe Ia identified 
without host galaxy spectroscopic redshift in DES that 
could be used for this type of analysis (Möller & the 
DES Collaboration in prep. 2024). 

'The DES supernova survey was supported by the 6- 
year OzDES survey on the Anglo-Australian Telescope 
(described in Lidman et al. 2020), which took multi- 
fibre observations of host galaxies to acquire redshifts 
of host galaxies of SNe. The total investment of this 
program was 100 nights, and for roughly 7596 of the 
targeted host galaxies a spectroscopic redshift has been 
secured. This program was fortuitous as the cameras 
for OzDES and DECam have a nearly identical field-of- 
view. Enormous resources would be needed to reproduce 
this joint program for LSST, which will find millions 
of SNe across 18,000 square degrees (Ivezić et al. 2019; 
Sánchez et al. 2022) (compared to the 27 square degrees 
of DES SNe). Surveys such as 4MOST will follow-up 
tens of thousands of these (Swann et al. 2019), but the 
full wealth of transient information may benefit from an 
entirely photometric approach. 

As statistical precision continues to improve thanks 
to the increased number of supernovae, a main theme 
for systematic analysis is second-order relations between 
different systematics. Typically, systematics are treated 
independently when building the covariance matrix. We 
have implemented a method to account for calibration 
systematics along with light-curve model systematics to- 
gether, but this is currently the only joint exercise. This 
type of work will grow in importance. For example, 
while photometric classification does not directly cause 
a large increase in the error budget, it hinders the abil- 
ity to constrain the intrinsic scatter model preferred by 
the data. Potentially, if LSST and other surveys such as 
those enabled by the Nancy Grace Roman Space Tele- 
scope have enough supernovae (Rose et al. 2021), the 
dataset can enable a forward modeling approach such 
as the Approximate Bayesian Computation method in- 
troduced in Jennings et al. (2016) and worked on in 
Armstrong et al. (in prep), which could vary all sys- 
tematics, nuisance, and cosmological parameters at the 
same time to compare against the data. 

Furthermore, as discussed in Section 5.1.2, modeling 
of the low-z sample remains a source of systematic un- 
certainty. This sample comes from a multitude of sur- 
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veys, even though we have removed many of the older 
inhomogeneous sources compared to analyses like Pan- 
theon+. In the near future, we expect additions from 
Zwicky Transient Factory (Smith et al. in prep. 2024), 
Young Supernova Experiment (Jones et al. 2021; Aleo 
et al. 2023), and Dark Energy Bedrock All-sky Super- 
nova Survey (DEBASS, PI: Brout) to improve low-z 
constraints of the SN Hubble Diagram, given their im- 
proved calibration and better understood selection func- 
tion. DEBASS will be particularly fruitful as it is a 
low-redshift sample taken with DECam, so a single in- 
strument and calibration catalog will be used for the 
full sample of DEBASS+DES, similar to the single- 
instrument PS1 sample in Jones et al. (2019). Using 
simulations, we estimate that quadrupling the size of 
our low-z sample (from ~ 200 to ~ 800 SNe expected 
from this next generation of low-z SN surveys) could en- 
able a reduction of uncertainties on w by ~ 30 per cent 
(for a FlatwCDM model, using SN data alone). 

Lastly, we note that while LSST and Roman may help 
improve a number of these issues, the first data release 
is still > 3 years away. We encourage work with the 
DES-SN sample as presented here, combined with other 
samples. Popovic et al. (2023b) recently showed the 
ability to combine separate photometric samples (PS1 
and SDSS) into the Amalgame sample (also shown in 
Fig. 12, and a similar analysis can be done by combining 
DES with these. It is reasonable to expect that with new 
low-redshift samples, and combination of high-redshift 
photometric samples, a sample with > 5000 likely SNe Ia 
can be compiled in the very near future. 


6. CONCLUSIONS 


The DES Supernova survey stands as a groundbreak- 
ing milestone in SN cosmology. With a single survey, 
we effectively tripled the number of observed SNe Ia at 
z > 0.2 and quintupled the number beyond z > 0.5. 
Here we present the unblinded cosmological results, and 
in companion papers make public the calibrated light 
curves and Hubble diagram from the full sample of DES 
Type Ia supernovae (Sánchez in prep. 2024; Vincenzi 
et al. 2024). 

After combining the 1635 DES SNe (of which 1499 
have a probability > 0.5 of being a SN Ia) with 194 ex- 
isting low-z SNe Ia, we present final cosmological results 
for four variants on ACDM cosmology, as summarized 
in Table 2. 


11 These upcoming low-z surveys are magnitude-limited rather than 
targeted, therefore they provide SN samples with a well defined 
selection function. 
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The standard Flat-ACDM cosmological model is a 
good fit to our data. When fitting DES-SN5YR alone 
and allowing for a time-varying dark energy we do see a 
slight preference for a dark energy equation of state that 
becomes greater (closer to zero) with time (wa < 0) but 
this is only at the ~ 2c level, and Bayesian Evidence 
ratios do not strongly prefer the Flat-wọwaCDM cos- 
mology. 

We compare cosmological results from each of our 
models to results from the CMB analysis of Planck Col- 
laboration (2020). There are some differences in the 
best fit values but in each case we find consistency to 
within 2c and a Suspiciousness statistic that indicates 
agreement among the datasets. 

Critically, the DES-SN5YR analysis shown here 
demonstrates that contamination due to SN classifica- 
tion and host-galaxy matching is not a limiting system- 
atic for SN cosmology; this opens the path for a new era 
of cosmological measurements using SN samples that are 
not limited by live spectroscopic follow-up of SNe. In- 
stead, our analysis shows the SN community that there 
are other factors that will be crucial for the success of 
future SN experiments: a high-quality low-redshift sam- 
ple, arobust UV and NIR extension of light-curve fitting 
models, excellent control of selection effects across the 
entire redshift range, and improvement in our under- 
standing of SN Ia intrinsic scatter properties and the 
role played by interstellar dust. 

Future work will conclude the Dark Energy Survey by 
combining these supernova results with the other three 
pillars of DES cosmology, namely baryon acoustic oscil- 
lations, galaxy clustering, and weak lensing. 


ACKNOWLEDGMENTS 


We acknowledge the following former collaborators, who have 
contributed directly to this work — Michael Childress, Ricardo 
Covarrubias, Chris D’Andrea, David Finley, John Fisher, Fran- 
cisco Förster, Daniel Goldstein, Santiago González-Gaitán, Ravi 
Gupta, Mario Hamuy, Steve Kuhlmann, James Lasker, Marisa 
March, John Marriner, Eric Morganson, Jennifer Mosher, Eliza- 
beth Swann, Rollin Thomas, and Rachel Wolf. 

T.M.D., A.C., R.C., S.H., acknowledge the support of an 
Australian Research Council Australian Laureate Fellowship 
FL180100168) funded by the Australian Government, and A.M. is 
supported by the ARC Discovery Early Career Researcher Award 
DECRA) project number DE230100055. M.S., H.Q., and J.L 
are supported by DOE grant DE-FOA-0002424 and NSF grant 
AST-2108094. R.K. is supported by DOE grant DE-SC0009924. 
M.V. was partly supported by NASA through the NASA Hub- 
ble Fellowship grant HST-HF2-51546.001-A awarded by the Space 
Telescope Science Institute, which is operated by the Association 
of Universities for Research in Astronomy, Incorporated, under 
NASA contract NAS5-26555. L.K. thanks the UKRI Future Lead- 
ers Fellowship for support through the grant MR/T01881X/1. 
L.G. acknowledges financial support from the Spanish Ministerio 
de Ciencia e Innovación (MCIN), the Agencia Estatal de Inves- 
tigación (AET) 10.13039/501100011033, and the European Social 
Fund (ESF) "Investing in your future" under the 2019 Ramón y 


Cajal program RYC2019-027683-I and the PID2020-115253GA- 
I00 HOSTFLOWS project, from Centro Superior de Investiga- 
ciones Científicas (CSIC) under the PIE project 20215AT016, and 
the program Unidad de Excelencia María de Maeztu CEX2020- 
001058-M, and from the Departament de Recerca i Universi- 
tats de la Generalitat de Catalunya through the 2021-SGR-01270 
grant. R.J.F. and D.S. were supported in part by NASA grant 
14-WPS14-0048. The UCSC team is supported in part by NASA 
grants NNG16PJ34G and NNG17PX03C issued through the Ro- 
man Science Investigation Teams Program; NSF grants AST- 
1518052 and AST-1815935; NASA through grant No. AR-14296 
from the Space Telescope Science Institute, which is operated by 
AURA, Inc., under NASA contract NAS 5-26555; the Gordon and 
Betty Moore Foundation; the Heising-Simons Foundation; and fel- 
lowships from the Alfred P. Sloan Foundation and the David and 
Lucile Packard Foundation to R.J.F. We acknowledge the Univer- 
sity of Chicago's Research Computing Center for their support of 
this work. 

Funding for the DES Projects has been provided by the U.S. 
Department of Energy, the U.S. National Science Foundation, 
the Ministry of Science and Education of Spain, the Science and 
Technology Facilities Council of the United Kingdom, the Higher 
Education Funding Council for England, the National Center 
for Supercomputing Applications at the University of Illinois at 
Urbana-Champaign, the Kavli Institute of Cosmological Physics 
at the University of Chicago, the Center for Cosmology and Astro- 
Particle Physics at the Ohio State University, the Mitchell In- 
stitute for Fundamental Physics and Astronomy at Texas A&M 
University, Financiadora de Estudos e Projetos, Fundação Car- 
los Chagas Filho de Amparo à Pesquisa do Estado do Rio de 
Janeiro, Conselho Nacional de Desenvolvimento Científico e Tec- 
nológico and the Ministério da Ciência, Tecnologia e Inovação, the 
Deutsche Forschungsgemeinschaft and the Collaborating Institu- 
tions in the Dark Energy Survey. 

'The Collaborating Institutions are Argonne National Labora- 
tory, the University of California at Santa Cruz, the University of 
Cambridge, Centro de Investigaciones Energéticas, Medioambien- 
tales y Tecnológicas-Madrid, the University of Chicago, Univer- 
sity College London, the DES-Brazil Consortium, the University 
of Edinburgh, the Eidgenóssische Technische Hochschule (ETH) 
Zürich, Fermi National Accelerator Laboratory, the University of 
Illinois at Urbana-Champaign, the Institut de Ciències de l'Espai 
(IEEC/CSIC), the Institut de Física d'Altes Energies, Lawrence 
Berkeley National Laboratory, the Ludwig-Maximilians Univer- 
sitat München and the associated Excellence Cluster Universe, the 
University of Michigan, NSF's NOIRLab, the University of Not- 
tingham, The Ohio State University, the University of Pennsyl- 
vania, the University of Portsmouth, SLAC National Accelerator 
Laboratory, Stanford University, the University of Sussex, Texas 
A&M University, and the OzDES Membership Consortium. 

Based in part on observations at Cerro Tololo Inter-American 
Observatory at NSF's NOIRLab (NOIRLab Prop. ID 2012B- 
0001; PI: J. Frieman), which is managed by the Association of 
Universities for Research in Astronomy (AURA) under a cooper- 
ative agreement with the National Science Foundation. Based in 
part on data acquired at the Anglo-Australian Telescope. We ac- 
knowledge the traditional custodians of the land on which the 
AAT stands, the Gamilaraay people, and pay our respects to 
elders past and present. Parts of this research were supported 
by the Australian Research Council, through project numbers 
CE110001020, FL180100168 and DE230100055. Based in part on 
observations obtained at the international Gemini Observatory, a 
program of NSF's NOIRLab, which is managed by the Associ- 
ation of Universities for Research in Astronomy (AURA) under 
a cooperative agreement with the National Science Foundation 
on behalf of the Gemini Observatory partnership: the National 
Science Foundation (United States), National Research Council 
(Canada), Agencia Nacional de Investigación y Desarrollo (Chile), 


Ministerio de Ciencia, Tecnología e Innovación (Argentina), Min- 
istério da Ciência, Tecnologia, Inovações e Comunicações (Brazil), 
and Korea Astronomy and Space Science Institute (Republic of 
Korea). This includes data from programs (GN-2015B-Q-10, GN- 
2016B-LP-10, GN-2017B-LP-10, GS-2013B-Q-45, GS-2015B-Q-7, 
GS-2016B-LP-10, GS-2016B-Q-41, and GS-2017B-LP-10; PI Fo- 
ley). Some of the data presented herein were obtained at Keck 
Observatory, which is a private 501(c)3 non-profit organization 
operated as a scientific partnership among the California Insti- 
tute of Technology, the University of California, and the National 
Aeronautics and Space Administration (PIs Foley, Kirshner, and 
Nugent). The Observatory was made possible by the generous 
financial support of the W. M. Keck Foundation. This paper in- 
cludes results based on data gathered with the 6.5 meter Magel- 
lan Telescopes located at Las Campanas Observatory, Chile (PI 
Foley), and the Southern African Large Telescope (SALT). The 
authors wish to recognize and acknowledge the very significant 
cultural role and reverence that the summit of Maunakea has al- 
ways had within the Native Hawaiian community. We are most 
fortunate to have the opportunity to conduct observations from 
this mountain. 

'The DES data management system is supported by the National 
Science Foundation under Grant Numbers AST-1138766 and 
AST-1536171. The DES participants from Spanish institutions 
are partially supported by MICINN under grants ESP2017-89838, 
PGC2018-094773, PGC2018-102021, SEV-2016-0588, SEV-2016- 
0597, and MDM-2015-0509, some of which include ERDF funds 
from the European Union. IFAE is partially funded by the 
CERCA program of the Generalitat de Catalunya. Research 
leading to these results has received funding from the European 
Research Council under the European Union's Seventh Frame- 
work Program (FP7/2007-2013) including ERC grant agreements 


19 


240672, 291329, and 306478. We acknowledge support from the 
Brazilian Instituto Nacional de Ciéncia e Tecnologia (INCT) do 
e-Universo (CNPq grant 465376/2014-2). 

'This research used resources of the National Energy Research 
Scientific Computing Center (NERSC), a U.S. Department of En- 
ergy Office of Science User Facility located at Lawrence Berkeley 
National Laboratory, operated under Contract No. DE-ACO2- 
05CH11231 using NERSC award HEP-ERCAP0023923. 


This manuscript has been authored by Fermi Research Alliance, 
LLC under Contract No. DE-ACO2-07CH11359 with the U.S. 
Department of Energy, Office of Science, Office of High Energy 
Physics. 


Facilities: CTIO:4m, AAT, Gemini:Gillett (GMOS- 
N), Gemini:South (GMOS-S), Keck:I (LRIS), Keck:II 
(DEIMOS), Magellan:Baade (IMACS), Magellan:Clay 
(LDSS3, MagE), SALT 


Software: numpy (Harris et al. 2020), astropy (As- 
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APPENDIX 


A. DATA RELEASE AND HOW TO USE THE DES-SN5YR DATA 


All the input/output files necessary to reproduce our analysis and the outputs of our analysis pipeline will be made 
available on acceptance of this paper. 

The DES-SN5YR analysis was run using the PIPPIN pipeline framework (Hinton & Brout 2020)!? that orchestrated 
SNANA codes for simulations, light curve fitting, BBC, and covariance matrix computation (SNANA, Kessler et al. 
2009),? and also integrated photometric classification from Möller & de Boissière (2020)! and Qu et al. (2021).1° 
Additional analyses codes that run outside the main pipeline include Scene Model Photometry (Brout et al. 2019b), fit 
to measure the SN population of stretch and color (Popovic et al. 2023a),1° SALT3 training (Kenworthy et al. 2021),!” 
and CosmoSIS to fit for cosmological parameters (Zuntz et al. 2015).!5 

We release the PIPPIN input files necessary to (i) generate and fit all the simulations used in the analysis (both the 
large ^biasCor" simulations to calculate bias corrections, and the DES-SN5YR-like simulated samples to validate the 
analysis); (ii) reproduce the full cosmological analysis, from light-curve fitting to photometric classification, distance 
estimates and cosmological fitting. Auxiliary files are also available within the SNANA library.!? 

The various (intermediate and final) outputs of our analysis pipeline will also be provided on github upon acceptance. 
This includes (i) light-curve fitted parameters, (ii) light-curve classification results, (iii) the final Hubble diagram and 
associated uncertainties covariance matrices, and (iv) the cosmology chains. 


12 https://github.com/dessn/Pippin 

13 https://github.com/RickKessler/SNANA 

14 https: //github.com/supernnova/SuperNNova 
15 https: //github.com/helenqu/scone 

16 https: //github.com/djbrout /dustdriver 

17 https://github.com/djones1040/SALTShaker 
18 https: //github.com/joezuntz/cosmosis 

19 https: //zenodo.org/records/4015325. 
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